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We propose a model for collisions between particles of a granular material and calculate the 
restitution coefficients for the normal and tangential motion as functions of the impact veloc- 
ity from considerations of dissipative viscoelastic collisions. Existing models of impact with 
dissipation as well as the classical Hertz impact theory are included in the present model as 
special cases. We find that the type of collision (smooth, reflecting or sticky) is determined by 
the impact velocity and by the surface properties of the colliding grains. We observe a rather 
nontrivial dependence of the tangential restitution coefficient on the impact velocity. 
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I. INTRODUCTION 

A rich variety of systems one encounters in nature may 
be considered as "granular gas" |lj . The most important 
difference between a "gas" of granular particles and a 
regular gas is the inelastic nature of the inter-particle 
collisions. The steady removal of kinetic energy from the 
granular gas due to dissipative collisions causes a variety 
of non-equilibrium processes, which have been subjects 
of experimental (e.g. E J E ffl [| ^ p[ p[ [TOf) and theo- 
retical (e.g. O, [l2|, [l3j [l4], pjjfl) interest. Particularly in 
recent time many of the experimental results have been 
reproduced and investigated using various techniques 
such as Cellular Automata (e.g. [Tt| , |l8| ), Monte- 
Carlo Methods JlStl, Lattice-Gas models [J20fl , Molecular 
Dynamics in two 




24 and three 
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dimensions and hybrid methods 

The loss of kinetic energy of a pair of inelastically col- 
liding grains can be described using the restitution coef- 
ficients for the normal and tangential components of the 
relative motion and e T 



-e N g N 
e T 9 T 



(0< 



< 1 



1< e T < 1) 



(1) 
(2) 



where g N and g T are the relative velocities of the par- 
ticles in normal and tangential direction before the colli- 
sion and (g N ) , (g T ) after the collision. 



Recently, the collision properties of small spheres have 
been investigated experimentally |52). These investiga- 
tions have shown that the type of the collision (sliding or 
sticking) depends on the ratio of g N and g T . The results 
were explained with different models for each type, and 
the coefficients in these models were independent from 
the velocity. On the other hand, laboratory experiments 
with ice balls |33| as well as with spheres of other materi- 
als (for an overview see J34j ) have shown that the normal 
restitution coefficient e N depends significantly on the im- 
pact velocity. As already seen, the tangential restitution 
coefficient depends on the impact parameters as well. 

The behaviour of the sheared granular material may 
be significantly different if the restitution coefficients de- 
pend on the impact velocity. This dependence should 
be taken into account in order to get an adequate model 
of the stress distribution |3q |. It is also known that the 
the parameters e N and e 1 crucially influence the global 
dynamics of granular systems (e.g. ]36], [37j). 

In the present study we investigate how the restitution 
coefficients depend on the relative impact velocity. For 
the normal component of the relative motion we derive 
an expression for the normal force acting between the 
colliding particles, which accounts for the dissipation in 
the bulk of material. One particular application of the 
results presented here is the explanation of experiments 
with ice balls 1 33 which are of importance for the in- 
vestigation of the dynamics of planetary rings B3]. A 



2 



static model for the inelastic impact of metal bodies was 
presented in |^| which based on the assumption of fully 
plastic indentation and constant mean contact pressure 



and led analytically to a proportionality e 



N 



[9 



-1/4 



for arbitrary material constants. In contrary, our qua- 
sistatic approach does not request other additional as- 
sumptions and can be adapted to different experimen- 
tal results by changing the coefficients in the differential 
equation which describes the time dependence of the de- 
formation. From these coefficients, material coefficients 
can be estimated . 

Our result contains the Hertz theory of elastic im- 
pact |4f[| and the theory of the viscoelastic impact by 
Pao |42fl as special cases. For the tangential component 
of the relative motion we consider a microscopic model of 
the contact of colliding particles. We derive a mean-field 
expression for the tangential inter-particle force. The re- 
sult contains the model of the tangential force of colliding 
particles by Haff and Werner jf3|, Q as a special case, 
and we are able to treat different tangential collisional 
behaviors within the framework of one single model. 

In section || we formulate the collision model and de- 
rive the equations for the normal and tangential compo- 
nents of the relative motion of the colliding grains. In 
section [II we present the results for the restitution co- 



efficients for the proposed model and discuss the depen- 
dence of the coefficients on the components of the impact 
velocity. A model for the dynamics of granular materi- 
als is proposed. In the last section IV we summarize 
the results. Details of the derivations are given in the 
appendixes |X| and [b]. 



one obtains Newtons equations for the translational and 
rotational motion 



dt 
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The force acting between the particles during collision 



consists of the normal component F^ = n 



and 



the tangential componenti^J = Fij — Fjj- Introducing 
the corresponding components g[- and of the relative 
velocity g^ and with 



_i _ 1 mjJi + m 3 Jj 

K ij — 1 H - - ~ 

./,./, : ;//, • ;// , • 
we rewrite eq. (0) omitting the indexes ij: 



f =F N /m et 
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(8) 



(9) 
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Using eqs. (|) the energy loss during the collision is 

AQ = ^ ((e N Y - (9 T ) 2 ((ST - 1 

The energy is conserved during the collision if e N = 1 
and e T = ±1. In these cases there is a completely elastic 
rebound for the normal component and either completely 
elastic rebound (rough spheres) or frictionless slipping 
(smooth spheres) for the tangential component. 



II. THE COLLISION MODEL 

We consider the inelastic collision between two spher- 
ical particles i and j. The values r^, Ri, fi, c3i, mi, and 
Ji are the position of the center of sphere i, its radius, 
velocity, angular velocity, mass and momentum of iner- 
tia, respectively. The relative velocity of the surfaces of 
the colliding particles at the point of contact is given by 
(e.g. p|) 



g*ij = [ri - uji x Ri rij - ifj + LUj x Rj rij 

= fi — fj — R4 uji x n — Rj ujj x n , (3) 

with n = — S 2 -, . Introducing the dimensionless moment 
of inertia Jj, the effective mass and the effective 
radius Rff 



Ji 



Ji 
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Reff = RjRj 



% ^ Ri + Rj 



(4) 
(5) 
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A Normal motion 

We assume that the colliding particles begin to touch 
each other at the time t = with the relative normal 
velocity g . When we introduce the deformation (or 
"compression" ) 



t(t)=Ri+Rj-(\m-m i) 

this velocity can be written as g N =| g N |= 
Thus from eq. (|^) we obtain the equations 
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N 
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m = g 

e(o) - o . 

The normal force F N consists of an elastic, conservative 
part due to the deformation £ of the particles and a vis- 
cous part due to dissipation of energy in the bulk of the 
particle material, depending on the deformation rate £. 
For the conservative part Hertz's theory of elastic con- 
tact H] gives for spherical particles 



(el) 
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where Y and v are the Young modulus and the Poisson 
ratio for the material the particles consist of. This rela- 
tion between the elastic component of the force and the 
deformation is valid for the quasistatic regime of the col- 
lision, i.e. when inertial and relaxation effects may be 
neglected (see appendix [b]). 

The existing phenomenological expressions for the dis- 
sipative part of the normal force which are either linear 
in the deformation rate £ (e.g. [[l3[ f45| ]) or quadratic fl6"| , 
however, do not agree satisfactorily with the experi- 
mental data for the normal restitution coefficient [ p3[ . 
Pao |Q extended the Hertz theory of impact for the 
viscoelastic case, where, however, the dependence of the 
bulk dissipation on the dilatation rate was neglected. In 
this theory memory effects in the dissipative processes 
were taken into account. Although the latter approach is 
not self-consistent (see appendix |B|), it predicts a power- 
law dependence of the dissipative force on the deforma- 
tion rate, yielding an exponent similar to that one for the 
quasistatic collision. In the present study we develop a 
self-consistent quasistatic approximation to calculate the 
normal force acting between colliding viscoelastic parti- 
cles. The quasistatic approximation is valid when the 
characteristic relative velocity of the granular particles is 
much less than the speed of sound in the material which 
is satisfied for many experimental situations even in as- 
trophysical systems such as planetary rings E7| . For the 
duration of the collision it is necessary to exceed sig- 
nificantly the viscous relaxation time in the material of 
colliding particles (see appendix |B|). 

Different from the approaches |34|, Q we take into ac- 
count both components of the dissipative force, arising 
from the shear strain rate as well as from the dilata- 
tion rate, which are both of comparable importance for 
the normal component of the relative motion. From the 
equation of motion for the viscoelastic continuum we find 
the general relation between the dissipative part of the 
normal force and the deformation rate. We show that 
memory effects in dissipative processes may be neglected 
in the case of a self-consistent quasistatic approximation. 
Since the calculation of the dissipative part of the normal 
force is rather straightforward, we present only the main 
idea of the derivation and refer to appendix [a| for further 
details. In appendix B the conditions for the validity of 
the quasistatic approach are given. 

The total normal force acting between viscoelastic par- 
ticles may be derived from a stress tensor combined of an 
elastic and a dissipative part |4q| 



with 



& = 0"(ei) + &(di 
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The displacements in the material is denoted by u and 
I is the unit tensor. -E1/2 and 771/2 are the elastic and the 
viscous constants of the particle material 



£1 = 



Y 



v 
Y 



3(1 - 2v) 



(18) 
(19) 



In the quasistatic regime (appendix |b|) the displacement 
field u (r, t) can be approximated by that of the static 
problem u(f). It is completely determined by the elastic 
component of the inter-particle force (|l4|). Thus, the 
displacement velocities can be written as 



• d 

u(r,t) ~ £ — u {el) , 



(20) 



where ur e n (f, £) is the solution of the static (elastic) con- 
tact problem. This expression depends parametrically on 
the deformation £ and the dissipative part of the stress 
tensor becomes 
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(21) 



The calculations can be significantly simplified when we 
notice that the elastic and the dissipative parts of the 
stres s ten sor are related in the quasistatic regime (see 
eqs. (WM and (^0|) 



>(dis) 



P 9 
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(E 1 <-> 77x , E 2 



(22) 



Therefore the impact problem for the viscoelastic parti- 
cles in the quasistatic regime can be mapped onto the 
corresponding problem for elastic particles. Performing 
calculations similar to that of the elastic case (for details 
see and appendix one can find an expression for 
the dissipative part of the normal force: 



A 



Y 



1 (37?2 - 771) 2 
3 (3772 + 277!) 



(23) 



Yv 2 



From eqs. (|23|) and ( J14| ) we obtain for the normal com- 
ponent of the relative motion 



3 77^(1^) (? /2+ 2 A ^t 



(24) 



with the initial conditions £ (0) = g N , £ (0) = 0. In the 
case of At; <§C £ equation (^4|) results from a Taylor- 
expansion of 



2 YVW' 



3 m e ff(l- v 2 ) 
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M) =0 



(25) 
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which formally coincides with the corresponding equation 
for the elastic problem, provided that £ is substituted by 

It has to be noted, that £ has its minimum at the 
beginning of the collision where £ takes its maximum. 
Hence, the condition A£ <C £ is not provided at the very 
beginning of the contact. On the other hand, the good 
confirmation of experimental facts (J33|) by the numerical 
solution of equation ( p5| ) points to its suitability for at 
least the rest of the collision time span. 

Taking into account (g N ) — £(i c ) (t c ~ the duration 
of the collision), the normal restitution coefficient is ob- 
tained from 



= i(t c )/i(o). 



(26) 



B Tangential motion 

In the idealized model the surface of contact between 
the spheres S is a perfectly flat circular area with radius 
Rs = y/2R e ff^(t). For the description of the tangential 
forces between the surfaces we follow a current model 
of tribology (e.g. |5(], |5~]] |) where the apparent surface of 
contact is built up of a large number of hierarchically 
ordered asperities varying in shape and size by several 
decades. For the processes of the momentum transmis- 
sion we will take into account only the largest-scale as- 
perities ("primary asperities"). The surface asperities do 
not affect the normal motion, if they are small enough 
(see appendix [A|), however they are responsible for the 
tangential forces, acting between the colliders. Here we 
consider a simplified mean-field approach and introduce 
the normal a N and tangential stress a T averaged over 
the contact area. Further we define the normal compo- 
nent of the total contact area of the asperities of both 
spheres S N which is responsible for the transmission of 
the normal force. Correspondingly the tangential projec- 
tion of the area S T is responsible for the transmission of 
the tangential force. These surfaces are related to the 
apparent contact area by the relations [p2 



S*( t )=f N (a N } S{t) 

s T (t) = f (>) S(t) , 



(27) 
(28) 



where the coefficients f N and f T depend on the average 
normal stress a N . When the spheres begin to touch each 
other, i.e. S = and = 0, we find f N (0) = and 
/ T (0) = 0. We expand the coefficients in eq. ( p8| ) with 
respect to a N = 0. The linear expansion yields for the 
tangential component of the surface 



s r (t) 



■ T a N S(t) 
df T 



(29) 



d<J N 



cr N =0 



For a given model of the sizes and shapes of the asperities 
one can calculate the value of 4> T |52] . In the case that 
the heights of the asperities obey a Gaussian probability 
distribution with mean value L one finds 



<j) T cx VI 



(30) 



For the average size of the asperities L of the surfaces the 
mean field approach yields the average shear deformation 
rj = b -j. The values C and b are the relative tangential 
shift of the particle-surfaces and a form-factor, respec- 
tively. We assume that the stress is uniformly distributed 
over the entire surface and find 



t t - 



-V = 



Y ( 



(1 + v) L 



(31) 



The linear relation between a T and r\ holds only for 
the elastic regime, i.e. only if a T does not exceed some 
critical value which is a specific material constant. If 
the shear stress exceeds this threshold <rj, the asperity 
which hinders the tangential relative motion of the sur- 
faces, is assumed to break resulting in a sudden release 
of the shear stress. At the same time the surfaces are 
shifted macroscopically with respect to each other by 



Co 



L (1 



b Y 



and one finds 



v (0 = v* 



(32) 



(33) 



n* 



Ko 
L 



where [^J denotes the integer of x. The breaking of 
the asperities dissipates the energy which was previously 
stored in the elastic stress, i.e. fracturing of the asperi- 
ties is the elementary dissipative process in the tangen- 
tial motion. From eq. ( |33| ) we obtain the shear stress as 
a function of the tangential displacement 



(0 
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) 



(34) 



and the tangential component of the inter- particle force 



F 1 





c 
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_Co_ 


) 
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c. 

Co 



(35) 



where F N — a N S is the normal component of the inter- 
particle force and /i = 4> T . 

It may be shown, that a more refined mean-field ap- 
proach, which does not use the assumption of the uni- 
formly distributed stress over the contact area, leads to 
the same eq. (K5J) for the tangential motion. 
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From eq. ([35|) follows the condition for the maximum 
tangential force 



F 1 



N 



(36) 



Thus our model reproduces the Coulomb friction law |53J 
with the friction coefficient /_t, expressed in terms of mi- 
croscopic parameters. The model for the tangential mo- 
tion is very similar to the extensively investigated one 
dimensional model by Burridge and Knopoff |55| in- 
tended to model earthquakes. 

i eff '((t) the 



With g T (t) = C(t), eq. (§) and F N = 



tangential motion is governed by the differential equation 



c-^w 



Co 







(37) 



with the initial conditions £ (0) 
value of £(t) is given by eqs. (|24| ) and (| 
tangential restitution coefficient reads 

e T = C(U/C(0). 



g T and C (0) = 0. The 
Then the 



(38) 



et al. 1984 - 
from eq. (18) 
from eq. (19) - 




normal impact velocity 



FIG. 1: The normal restitution coefficient e N versus the 
normal component of the impact velocity g N measured in 
cm s -1 according to eqs. (18) and (19). The dashed line 
denotes the dependence e N (g N ) measured by Bridges et 
al. [33] 



III. RESULTS AND DISCUSSION 



The obtained equations for the normal (eq. (|25|), (j2J)) 
and tangential motion (eq. ((3^)) have been solved nu- 
merically using a Runge-Kutta Method of fourth order 
with adaptive step size J5(|. The restitution coefficients 
e T and have been calculated as functions of the nor- 
mal and tangential relative velocities g T , g N . For the 
integration we used the parameters of ice at low temper- 
atures [[57]]: Young modulus Y — 10 GPa, the Poisson 
ratio v = 0.3, the particle size R — 10~ 2 m with density 
p = 10 3 kg m~ 3 . The coefficient A in eq. ( pi| ) was con- 
sidered to be a fit parameter, due to lack of information 
about the dissipative coefficients 771 and 772 . Fig. |l| shows 
the numerical result of our model for the normal resti- 
tution coefficient e N as a function of the normal relative 
velocity g N compared to experimental data for the colli- 
sion of spherical ice particles with an ice wall |53|. The 
experimental results are well reproduced by our model. 

For the investigation of the tangential restitution co- 
efficient of colliding homogeneous spheres (J = |mi? 2 , 
k = |) we have chosen the Coulomb friction coefficient 
from the interval /j, € [10 -2 ... 11. The value of erj is a 



material constant. With eq. (30) and the definitions of £0 
and 77, we estimate Co which characterizes the size of the 
surface asperities via /i = a -\/Co- In the numerical cal- 
culation we have chosen a — I. The results are shown in 
fig. ||. The tangential restitution coefficient e T is drawn 
versus the plane defined by the tangential and normal ve- 
locities g T and g . The three plots correspond to the val- 
ues of the asperity sizes Co = (l0~ 7 ; 2 • 10~ 4 ; 10~ 3 ) R e ff, 
respectively. 

The obvious common feature of all cases is sliding of 
the surfaces (e T > 0) for small g N and large g T . This 



is quite plausible since smaller impact velocity g N corre- 
sponds to a smaller normal acceleration and thus, to a 
smaller value of the maximal tangential force, eq. d35|). 
As a result e T — > 1 at g N — > due to vanishing tangential 
acceleration. At the same time, for the high tangential 
velocity, g T (0) ^S> 1, (g N (0) ~ 1) sliding occurs owing to 
a considerable breaking of the asperities. The area of the 
sliding phase in the g N -g T plane depends on the size Co 
of the asperities. 

In the case of Co = 10~ 7 R e $ sliding occurs in the entire 
velocity range according to values 0.85 < e T < 1. The 
small asperities are not able to cause a sufficient torque to 
change considerably the spin of the individual particles. 
Here we are close to the case of ideal smooth spheres 
where no change of the tangential motion is expected 
(e T = 1). 

In the other two cases Co = (2 • 10~ 4 ; 10~ 3 ) R eff one 
recognizes two phases: 

1. Sliding e T > at small g N and high g T 

2. Reversal of the spin of either particles e T < at 
small g T and higher g N . 

Case |] corresponds to the effect discussed in the context 
of Co = 10~ 7 . Despite of being far from rather smooth 
spheres, the small tangential force originated from small 
g N , changes the velocity g T only slightly. Hence one 
has e T > 0, which is also the case for high velocities 
g T where the asperities break. In the case || we have 
the other extreme: a high normal acceleration causes a 
tangential force which is high enough to change the sign 
of g T as long as the asperities do not break (small g T ). 
A complete reversal of the tangential velocity according 
to g T (0) — > — g T (t c ) is not possible because of the 
dissipation arising of the bulk viscosity of the material 
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FIG. 2: The stereographic projection of the tangential 
restitution coefficient e T versus the plane g N -g T of the 
tangential and normal components of the impact velocity. 
The three parts of the figure belong to different values 
of the "size" of the asperities: (a) Co = 1CP 7 R e ", (b) 
Co = 2 • I(T 4 R e ff, (c) Co = 10~ 3 R eff . 



which enters the normal as well as the tangential forces 



(see eqs. ©,©). 

Both types of behavior of the tangential motion are 
separated by a sharp transition at g T — gj cr ^ where 
the asperities begin to break (see surface plots for Co = 
(2-ICT 4 ; ltT 3 ) R e B). The higher g N the larger the 
critical tangential velocity gT cr \ ■ A higher normal veloc- 
ity g N causes a stronger counteracting force F T and thus 
a larger tangential impact speed g T is necessary to reach 
the critical deformation where the asperities break. Both 
cases (Co = (2 • I0~ 4 ; 10~ 3 ) i? e -"j reveal similar qualita- 
tive behavior but the ranges of different type of motion 
(0, cover different areas in the g T -g N plane. 

The results show that our model includes a continuous 
transition from the limit case of rough spheres (e T — ► — 
1) to the limit case of smooth spheres (e T — * I). In 
the literature of the dynamics of granular material an 
alternative step function is widely used for the tangential 
force E§ 



-j s m " \g 



(39) 



The numerical evaluation of the considered model 
(fig. |2|) reveals surprising behavior of the tangential resti- 
tution coefficient e T as a function of the normal velocity 
g N at fixed tangential velocity g T . (This effect is no- 
ticeable for the largest values of Co-) At low and mod- 
erate g N , e T first decreases with increasing g N down to 
its minimal negative value in a manner discussed above, 
but when g N exceeds some threshold (approximately of 
several g T ), it starts to increase up to zero at very high 
values of g N . This effect may be explained as follows: 
For high values of g N the average normal force is large 
and causes thus a large tangential force, which can effec- 
tively decelerate the initial tangential velocity without 
switching to the sliding regime. 



Calculating the restitution coefficients e 



(in the 



limits of our model) we obtain a complete description 
of binary collisions. Therefore one can determine the 
dynamics for moderately dense granular gases, where an 
evolution occurs via a sequence of binary collisions. For 
such systems we have the following Boltzmann equation 
for the one-particle distribution function 



^ + $1 ■ Vj /(I) = J dv 2 J dto 2 J dn \g N ■ n\ 9 



9 N • n) 



/(1/)/(2 2 -/(l)/(2) 



(40) 



with 0(a;) given by 

and with the common notations, e.g. (1) = {ri,Vi,Ui,t}. 



The velocity and angular velocity of the first particle after 
the collision v[ and u[ can be expressed in terms of the 
pre-collisional values via the relations 



< = € + ^ { V T ( 9 \ 9 T ) ~ 1] 9 T - [e N ( 9 N , 9 T ) + 1] 9 N } 



(42) 
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i<0 



0J-, 



..0 1 + g r RHx{ [? {g\g T ) - 1] f - [e» (g\g T ) + l] g N } 



(43) 



and analogously for v^, lu^- With the use of eqs. ( f4(i| ) 
and (H) and the above calculated restitution coefficients 



e " (3" >9") an< i e T (g N 1 9 T ) ( e Q s - (26) and J3S|)) one can 
describe the evolution of moderate dense granular gases 
without computing the detailed dynamics of binary col- 
lisions as it is usually done in the "soft sphere" Molec- 
ular Dynamics (MD) approach. Here one considers the 
grains as elastic bodies which deform each other during 
a collision. There are several ansatzes for the force act- 
ing between touching grains In all cases 
one has to chose a time step for the integration scheme 
which is significantly smaller than the typical collision 
time. Hence, during each collision one has to calculate 
about 10. . . 1000 times the interaction force between the 
grains to provide satisfying accuracy of the simulation. 
When two grains approach each other they do not feel 
any interaction as long as they do not touch each other. 
When granular particles collide they interact via huge 
restoring forces which can be expressed by Young mod- 
uli of the order of Y = 10 7 kg/m sec 2 . The difficulty 
of the simulation consists in the extreme short range in- 
teraction of the particles and the resulting huge gradient 
of the interaction force. Therefore presently one cannot 
simulate much more than 3000 granular particles in three 
dimensions (e.g. [5£, |6(|) and about 10 4 particles in two 
dimensions (e.g. [31]). 

Another method for the simulation of granular assem- 
blies is the "hard sphere" approach where one does not 
consider the details of the collision but only the pre- and 
post-collisional velocities of each pair of colliding grains. 
The goal of these simulations is the low numerical com- 
plexity. One needs only computational effort when parti- 
cles collide but not in between the collisions. This allows 
for the application of so called event-driven calculations 
(e.g. Ip^l ). Hence, one can simulate much more particles 
than with "soft particle" methods. 

One of the preconditions for the application of the 
"hard sphere approach" is the exact knowledge of the 
normal and tangential restitution coefficients, e N and e T , 
as functions of the normal and tangential impact rates, 
g N and g r , which theoretically to determine was the goal 
of the present paper. 

An interesting possible application of this approach is 
the dynamics of planetary rings composed of icy and sili- 
cate material which is determined by inelastic dissipative 
collisions . The calculation of such systems using the 
traditional MD is impossible due to the huge number of 
particles in these systems. 



IV. CONCLUSION 

A novel model for collision of particles in granular gases 
is proposed. For the normal component of the relative 



motion the equation of motion is derived based on the 
general consideration of the viscoelastic impact. We find 
the expression for the dissipative part of the normal force 
in the self-consistent quasistatic approximation which 
generalizes the existing results for the viscoelastic col- 
lisions f42fl . For the tangential relative motion we inves- 
tigated a microscopic model of surfaces of the colliding 
particles which are in contact. We found a mean-field ex- 
pression for the tangential inter-particle force, which can 
reproduce smooth, reflecting or sticky collisions depend- 
ing on the microscopic parameters of the surfaces and on 
the relative impact velocity. A frequently used model for 
collisions of granular particles by Haff and Werner 
is contained in our model as a special case. The restitu- 
tion coefficients for the normal and tangential motion are 
calculated as functions of the relative impact velocity. A 
rather nontrivial strongly non-linear dependence of the 
tangential restitution coefficient on the impact velocity 
is observed. 

The obtained restitution coefficients may be used to 
describe the dynamics of moderately dense granular 
gases, where the evolution occurs via a sequence of suc- 
cessive binary collisions. 
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APPENDIX A: GENERALIZATION OF THE 
HERTZ THEORY 

We briefly sketch the Hertz theory of elastic impact 
and give a generalization of this theory for the case of 
viscoelastic collisions. 

In the quasistatic approximation which is used in 
Hertz's impact theory it is assumed that the (time de- 
pendent) strain and the (time dependent) stress are re- 
lated in the same manner as in the static case. It may 
be shown (see appendix [b]) that this approximation is 
valid for the elastic case when the characteristic velocity 
is much less than the speed of sound in the material of 
the colliding particles. Moreover for the viscoelastic case 
it is required that the viscous relaxation time of the ma- 
terial is much shorter than the duration of the collision. 
In the static case the equation of equilibrium reads [M 







(Al) 



where the elastic stress tensor <7 ( e i\ i s expressed in terms 



of displacements uir) via eq. (HA). Hence the static 
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eq. (Al) can be written as 



V 2 ± u + b 2 Vjju = 
fo2 _ AEi + 6E 2 _ 2 (1 - v) 



(A2) 



3Si 



(1-21/) 



with the "longitudinal" and "transversal" parts of the 
Laplacian. 



Vj=VoV 

v 2 = v 2 - v 2 



(A3) 
(A4) 



Th e b oundary conditions for the displacements in 
eq. ( [A2] ) are formulated on the surface of contact. From 
geometric considerations it follows that the contact area 
between two colliding particles is a plane. Using the ap- 
propriate coordinate system centered in the middle of the 
contact region (where we set z = 0) one can write 



Cix 2 + C 2 y 2 + u z \ + u z2 = i 



(A5) 



The values u z i = u z x(x,y) and u z2 = u z2 (x,y) are the 
^-components of the displacements in the materials of 
the bodies at the plane z = 0, £ is the total deformation 
(the sum of the deformations of both bodies at the center 
of the contact area, i.e. at x = y = 0). The constants C\ 
and C 2 are expressed in terms of radii of curvature of the 
surfaces in contact (see e.g. JO], [l8|). The values of u z x 
and u z2 may be expressed in terms of the normal pressure 
P z [x,y) that acts between the bodies at z — P8fl : 



A 

u zl {x,y,0) = u z2 = — 

7T 



dx'dy' (A6) 



A 



= \J{x - x 1 ) 2 + (y - y'Y 
2Ei + 3E 2 1 - 



E 1 (Ex + 6E 2 ) 



Y 



For simplicity we assume that the colliding particles are 
of the same material. The normal pressure P z is related 
to the total normal force Fi 



(el) 



3 Ft, 



P{x v) - - - («o ,/i x2 y 2 

z[X,V) 2irab V a 2 6 2 



(A7) 



where a and b are the semi-axes of the contact ellipse. 
The latter values as well as the compression £ may be 
found from the set of equations 



(e0 



-A 



dq 



n 2 Jo vV 2 + v) (b 2 + q)q 

dq 



r- o />oc 



c 2 = 



tt 2 J (a 2 + q) ^(a 2 + q) [b 2 + q) q 



(A8) 
(A9) 



7(ei) -^A 



dq 



7T 2 J (& 2 + g) ^(a 2 + 9) (b 2 + q) q 



(A10) 



From the above expressions follows that for all bodies 
in contact having smooth surfaces (in the mathematical 



sense) the total force and the deformation are related via 
the power law 



Fun (0 = c e' 2 



(AH) 



The constant c depends on the elastic properties of the 
materials and on the local curvatures of the colliding bod- 
ies. For the case of the spherical particles one has the 
Hertz's law 



F(el) (0 



2 Y 



3 



(A12) 



Using this relation between force and deformation and 
the equation of motion (eq. (|l5|)) one can describe the 
elastic collision completely. The duration of the collision 

is ([0 B) 



ic = 2 <M f^f V)) 



4 2 



(A13) 



e= . 
V53A/ 

In the solution of the elastic contact problem the dis- 
placement fields u\ (r) and u 2 (r) are completely defined 
by the value of F( e/ ) and thus by the value of the de- 
formation £. Hence we write u(r) = u(r, £), i.e. the 
displacement field depends explicitly on the compression. 
Therefore we obtain for the velocity of the displacement 
in the quasistatic approximation 



u (r) 



C^(f,0 



(A14) 



and correspondingly for the dissipative part of the stress 
tensor 



'(dis) 



> d_ 



r\\u ik 



{Ex ^m, E 2 ^ m ) (A15) 



We emphasize that the expression in the curled brackets 
in the above equation coincides with the elastic stress, 
provided the viscous constants are substituted by the 
clastic ones. The latter expression for the dissipative 
stress tensor is written for the case when the memory ef- 
fects in the viscous processes may be neglected. A more 
general case is discussed in appendix |B]. 

The a^ n component of the elastic stress is equal to the 
normal pressure P z at the plane z = 0, 



<ei)(^y:°) 



du z 
oz 



3 F 



(el) 



2 nab 





( du x 


=L _ 


du z 


-f) 


\ dx 


dy 


V 

dz 



1 - 



b 2 



(A16) 
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With the transformation of the coordinate axes 



V = oty 
i 



and 



V2 + lvJ \E 2 ~^E 1 



a (E 2 - §£i) 



we obtain 

du z 



dz 



a = aa 
b = ab' 



( 771 \ / du x du y du 



(A17) 
(A18) 
(A19) 



(A20) 

(A21) 

(A22) 
(A23) 



a e 



duz 
dz 1 



E- 



dx dy dz 
_ dux_ duy_ du z 

3 J\ dx' dy' dz' 



.3F, 



(el) 



x' 2 y' 2 
2 na'b' V ~ a/* ~ V 2 



2 7ra6 



x 2 y 2 
a 2 6 2 



(A24) 



Applying the operator £ ^ to the previous expression we 
obtain the result for the viscous stress. Integrating the 
viscous stress over the contact area we finally find for the 
dissipative component of the inter-particles force 



• d 

E(dis) = AS.-^F^ (£) 



(A25) 



A = a?0 = i;;?7^ {1 -*1 { }- 2V) -(A26) 



3(3772 + 277!) 



Yv 2 



Thus one obtains for the normal force which acts be- 
tween the viscoelastic bodies in the quasistatic regime of 
collision 



F = const f e /2 + i 



(A27) 



The constant in eq. (A27) coincides with that for the 
elastic force. For colliding spherical particles we arrive 
at eq. (||). 

The impact theory, sketched above, was developed for 
the bodies with the smooth surfaces. If the surface asper- 
ities are taken into account, one can consider the actual 



surface as a smooth one (obtained by averaging over the 
asperities heights), with a small perturbation superim- 
posed due to the presence of the asperities. One can 
also consider the actual normal displacements and nor- 
mal pressure as a sum of the averaged (over the asperi- 
ties) values and the small perturbation. Then it is easy to 
show that the equations, obtained for the averaged values 
coincide (due to linearity of the problem) with the corre- 
sponding equations for the elastic collision of the smooth 
bodies. As a result, the relation between the force and 
deformation, £ is the same as in the Hertz's theory, pro- 
vided that £ is defined with a use of the averaged over 
the asperities radii of the colliders. 

Considering the normal motion for the dissipative colli- 
sions, one need not to consider the plastic deformation of 
the asperities, since the size of the asperities is assumed 
to be very small compared with the radii of the spheres. 
For our calculations in fig. 2. the asperity size is 10 3 to 
10 7 times smaller than the effective Radius of the parti- 
cles. Hence the dissipation in the bulk of the asperities 
is negligibly small compared to the total dissipation in 
the compressed part of the collider. Moreover, the ratio 
of the normal to tangential stress, may be roughly esti- 
mated as, a N ja T ~ (^/i?) 1 ^ 2 , so that the crushing of the 
asperities does not seem to be important for the normal 
motion, if £/i2 1 and if the conditions of the qua- 
sistatic collision hold. Thus one concludes, that the sur- 
face asperities may be ignored, when the normal motion 
is studied, provided they are small and the conditions of 
the quasistatic collision are satisfied. 



APPENDIX B: VALIDITY OF THE 
QUASISTATIC APPROXIMATION 

To analyze more rigorously the conditions when the 
quasistatic approximation is valid we write the equation 
of motion for the viscoelastic continuum 



pu = V 



>-( e i) 



>{dis) 



(Bl) 



where p is the density of the material. The expressi on fo r 
the elastic part of the stress tensor is given by eq. ( |ll A| ) . 
Taking into account the memory effects of the dissipative 
processes in the material one can write for the dissipative 
part 



V(dis) (t) = Ei I tpi(t- t) 



J ipi(t-r) ^ jv o U{t) + ff(r) o v} - ^ IV • U{t) + E 2 tp 2 (t- r) /V 



u(r) (B2) 



where the (dimensionless) functions an d fait) are coincides with the corresponding expression for the vis- 

relaxation (or "memory") functions for the distortion cous stress tensor in |34|, for tp 2 (t) — 0. The latter 
strain and the dilatation, respectively. Note that eq. (B2) approximation means that one neglects the bulk dissipa- 
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tion due to the dilatation rate. For the normal motion 
of colliding particles, however, the dissipation of energy 
due to the dilatation rate and the dissipation due to the 
distortion strain rate are of the same order of magnitude. 
Thus we keep both relaxation functions in our considera- 
tions. Introducing transversal and longitudinal velocities 
of sound in the material 



Y 



c t = 



Eh. 



2p 2p{l + v) 



(B3) 



2E, + 3E 2 



Y(l-v) 



3p p (1 + v) (1 - 2v) 

2(1-1/) 



(B4) 



(1 - 2v) 



one can write the equation of motion for the viscoelastic 
medium 



1 



' Vjjtt 



| + viVi vjj 



-ipi * u + b 



(B5) 



where ip * u denotes convolution. 

T o es timate the relative importance of the terms in 
eq. ( p35| ) we introduce the characteristic velocity vq = 
g N and the characteristic time tq = t c , where t c is the 



duration of the collision, introduced above in eq 
Then the characteristic length is Rq = vo tq. Eq 
can be written in a dimensionless form 




^ t(t) = {Vi«(t) + 6 2 Vfiu(f)} 



T vis,l \ r-i2 ~ 

V i U 



t vis ,i/2 = / 1P1/2 (T)dr . 
Jo 



(B7) 



We use the following representation of the convolution 



if;*u = u(t*) / Vi/2 r ^ T 
J o 



(B8) 



is of the order of 1, while by the definition of t. 



'vis, 1/2 



r vis, 1/2 



ipi/2 (r) dr 



and consequently 



1p*U~ u(t)T viSj i/2 



(B9) 



(BIO) 



1 > 4 

Ct 



1 > 



,1/2 



70 



(Bll) 
(B12) 



with in being a dimensionless time from the inter val < 
<* < t = t/rQ. The relation for the convolution (B8) is 
valid if ip(t) > 0. 

During the collision process t is of the order of tq, i.e. t 



these 



values are of the order of the relaxation times for the 
dissipative processes in the material. That means, that 
T v is, 1/2 characterizes the time when the memory effects 
are important. If the duration of the collision is much 
greater than the relaxation times, i.e. if T vis ^/ 2 <C t , 
one can write 



If the characteristic velocity vq is much less then the 
speed of sound in the material too, one can neglect the 

terms with vanishing factors («o/ c t) an( i i n 

cq. (B6) and finally one arrives at the static eq. flA2p. 



That means the quasistatic approach is valid provided 
that the conditions 



7~vis,l 
TO 



4 ^ 
Vjj < -u (U 



b 2 -- 



(**) 



(B6) 

hold. From the above considerations follows that in the 
quasistatic approximation the memory effects in the dis- 
sipative processes are not important and the viscous part 
of t he str ess tensor may be written in the same way as in 
eq. (HA), with the viscous constants rji and 772 given by 



V1/2 = E1/2 T viS; x/2 = E1/2 I ipi/2 (t) dr (B13) 

Jq 

It is worth to note that the quasistatic approximation 
is valid for many of the granular gases one encounters in 
nature, since usually the characteristic velocity in these 
systems is low. One should also note that the description 
of the collision in the quasistatic approximation is rigor- 
ous in a sense that no other additional approximations 
are used. 

As it follows from the above considerations, it is not 
correct to use the time dependent relaxation functions 
for the dissipative part of the stress tensor together with 
Hertz's quasistatic relations since this approach 

is not self-consistent and one has to assume a lot of ad- 
ditional hardly controllable approximations. 
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